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1 Introduction 



In the following we will be concerned with the calculation of massive loop inte- 
grals. We will present methods to calculate such integrals for the 1- and 2-loop 
case. For the 1-loop case these methods allow for the calculation of arbitrary rank 
tensor integrals for arbitrary n-point functions with arbitrary integer powers of the 
propagators. The methods for the 2-loop case cover the scalar planar and nonpla- 
nar integrals. The tensor cases can be handled by the usual reduction to 1-loop 
integrals or, in a much more direct approach, by integrating the 'characteristic 
integrand' of the graph under consideration 

In this article we give only typical introductory examples of these methods. We 
refer the reader to the literature to study the methods in detail || |o[ 0, [H], |l2[ . 

2 1-loop Integrals 

2.1 The Scalar Function 

In this section we calculate the two-point function 

1{q2) = I d ° l ik' (1) 

where 

Pi = I 2 - m\ + ig 

Pi = {1 + q) 2 - m 2 2 + ig. 

Here we do not specify the masses mi, TO2, so we want a result which is valid for the 
massive as well as for the massless case. We work within dimensional regularization 
in the sense of Q, D = 4 — 2s is the complex dimension. 

There is only one scalar-product involved in the integrand Eq.(|l|) and so if we 
choose to be in the rest frame of q = qo, we can separate the Zo-hrtegration. No 
non-trivial angular integration is then involved in the (D — l)-dimensional space- 
like integration. So we can do D — 2 angular integrations trivially, leaving the Iq 
and l± integrations to be done. We have denoted the modulus of the space-like 
part of I by l±. 



The integral Eq.(|l]) then becomes 
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where V\ = l — l ± — m\ + ig 

Vi = {l + go) 2 - lj_ - m\ + ig. 
The ^-integration in Eq.(|^) can now easily be performed. The result is 

1 



I{q 2 ) = C dl 



,l-2e „..l-2el 



where w\ = J —Iq + m\ — ig 



u> 2 = 



\J-{Iq + q) 2 +m\-ig 



D-1 

7T 



C = T ^_r(l/2 + e). (3) 



1 



Now the remaining /o-integration becomes 



Hi 2 ) 



M d = 



Ce {-i*{l/2-e)) 

i - ml 



dl 



(I 2 - m\ + ig)* 

q/2 + M d + l 



q/2 -M d -l 



2q 



(4) 



where we used translation invariance for the second term in Eq. ([}]), translation 
invariance being guaranteed by the very definition of dimensional regularization 
(Re(e) big enough in intermediate steps). The integral Eq.(^) can be determined 
by various methods. We will do it in the following by introducing 7?.-functions 
H 11 which are the typical ingredient of our 1-loop method. 
We modify the integral Eq.(|]) as follows. 



I = 



Ci exp(«7re) 
2^ 



OO J 

as 



(q/2 + M d ) 



(s — m\ + ig) 2 2e 
-{q/2 + M d y + s 



+(q/2 - M d ) 



(s — m\ + ig) 
-(q/2-M d y 



i-2e 



(5) 



which gives in terms of 7?.-functions 
i exp(i7re) C 



I 



2q 



(q/2 + M d )B(-,e) 
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+ e. 1; ~™\ + iQ, -(9/ 2 + M df) + {q/2 - M d )x 
xB(i e)n. e {~ + e, 1; -m\ + ig, -(q/2 - M d ) 2 ) 



(6) 



The Beta- function in Eq.(g) involves a pole in e. Expanding the following TZ- 
functions to 0(e) (see ||) we can express our result in terms of logarithms. The 
result is 



I{1 2 ) 



- - 7 - logTT + 2 + -(1 + m\/q 2 - m 2 /q 2 ) x 



1 2 

z\ log — - — ; logmf 



1 + Zl 



+ ^~rn 2 /q 2 + m 2 /q 2 ) X 



where 



1 — (1 — z 2) , 01 

3a log — — logm|J 

1 + z 2 



Zl = a 1 - 



Am\q 2 



(q 2 + m 2 — m 2 .) 2 

Arn^q 2 
(q 2 — m 2 + in 2 ) 2 



Z2 



where the upper (lower) sign has to be used if the argument is in the right (left) 
complex half-plane corresponding to Carlson's conventions on the scaling law for 



2 



7?.-functions and quadratic transformations on them j3j. This is the well-known 
result which one also obtains by conventional methods. 

2.2 The Tensor Functions 

As it is shown in || the above reduction allows also to calculate arbitrary tensor 
integrals. Here we give only the result. 

I ]k {q 2 ) = C(B(e- k/2 - [j}/2, 1/2 + [j]/2)(q/2 + M d y+^ 
x (-l)^ft_ e+fe/2+[j]/2 (-l/2 - k/2 + e, 1; z uyi ) 



i=0 



( J ) (-qy- l B(e - k/2 - [i]/2, 1/2 + [i]/2)(g/2 - M d f+ 1 ^ 



where 



x^_ £+fe/ 2 + [ s :]/ 2 (-l/2 - k/2 + s, 1; z 2 , 2/2 ) ) 
rn, to even 
to + 1, to odd 

^ 7r"^2^" iexp(i7re) (— l) fc / 2 

' _ l/2 + fc/2-e 2^ X 

T(3/2 - £ + fc/2)r(l/2 - k/2 + e) 
X r(3/2-e) 
Z\ := -TOi + yi = -(q/2 + M d ) 2 
z 2 := -to| + ig,y2 = -(q/2 - M d ) 2 . 

In the above formula for Ij k (q 2 ) j denotes the number of indices in the parallel 
space and 2k the number of indices in orthogonal space. One then has for a tensor 
integral 

I^... H = 0, j 2 odd 

= Ij 1 j 2 N (h)T Hl+1 ... H , h even, 

where 

T Hl+1 ... H = g^ lHi+2 ■ ■ ■ gf?~l H + a11 permutations of {fi jl+i . . . flj} 

N(j 2 ) = g H ' +lH ^ ■■■g H - lH T Hi+1 ... H . 

Here we assumed without loss of generality that the first j\ indices lie in the parallel 
space. We omitted these indices on the right hand side of the above equation 
because they all have fixed value fj,i = ... = fj,^ = (the parallel space is one- 
dimensional in this simple example). So, the expression on the right hand side is a 
j 2 -rank tensor in the orthogonal space. This tensor structure is given by the totally 
symmetric tensor T, constructed as the symmetric product of metric tensors in the 
orthogonal space. N(j 2 ) normalizes the tensor T. With 1 we denote this metric 
tensor in the orthogonal space. One has g®^ 1 ^ — g^T^g^" = D — 1. 

In the derivation of this formula the celebrated 'existence theorem of associated 
functions' has been used ||. 

The calculation of higher powers of propagators can now be easily done by 
applying appropriate differential operators with respect to masses to the result 
expressed in 1Z- functions || . 
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The next step in our examination of one-loop integrals is the calculation of the 
three-point functions. These form factor functions have infrared and on-shell singu- 
larities. They are, therefore, a more sophisticated test of our method. Nevertheless 
they can be calculated following the same route as for the 2-point functions flC| ]. 
The same is true for arbitrary n- point functions I7L 13 1. 



In the following we shortly recapitulate the calculation of the 3-point functions. 
We use the following notation for the one-loop three-point scalar integral 

'= /<*%HUr. (7) 



P0P1P2 



where 



Pa = 11-11-11-1^1 + ir] 

Pi = (lo + qwf -lf-l±-mf + iri 

P2 = (h + <?2o) 2 - (h + <?2i) 2 - l± - mj + ir), 

so that Iq and l\ span the parallel space (the now two-dimensional linear span of 
the exterior momenta), while l±, as before, is the modulus of the loop-momentum 
in the orthogonal space. 

The resulting non-trivial integrations are 

dl / dh / dl ± -±—, (8) 

-00 J -00 JO fQf\fl 

where C= 2^ /T(^). 

With the help of a partial fraction we have three summands to calculate. The 
^-integration is easy to do. The next integration can be done via the residue 
theorem and the remaining integrals will again be interpreted as 1Z- functions. 

Introducing the abbreviations 



aoi 


= -2gio 


601 


= 


coi 


2 2 

= m 1 — m — 


?io 


002 


= -2(720 


b{>2 


= 2g 2 i 


C02 


2 2 
= m 2 -m - 


gio + gli 


«ii 


= 2gi 


hi 


= 


Cll 


= -ml + ml 


- g?o 


ai2 


= 2(gio - g2o) 


l>12 


= 2g 2 i 


C12 


= —m\ + m\ 


- (gio + g2o) 2 + gli 

- gio + gli 


021 


= 2g 2 o 


l>21 


= -2921 


C21 


= -ml + ml 


a22 


= 2(g 2 o - gio) 


1)22 


= -2g 2 i 


C22 


= -ml + m{ 


- (gio + g2o) 2 + gli, 



wc find the final result to be 



2 4 A -e 

I = CB(2 £ ,lW-£,e)EE a '# ( 9 ) 

i=0 1=1 l ' 1 



TZ-2e(£, e, 1; -vti, -Vii , A,i/Ci,0 



where 

vti 



Bid 




mf — irj 




V 4A h 


Ai t i 


Bid 




mf — ir/ 


2A iA 


V 4 4< 


Ai,i 



Oil = 



-1,1 = 1,2 
-1,1= 3,4 
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aj,i + frj.i 
o-i,i — bi,i 

Qi,2 + bj,2 

Cj,l 

Cj,2 
Cj,2 — &i,2 

ai,i(ai,2 - fy.a) + a*,2(6i,i - a»,i) 

Cti,2(ai,l - + Cti,l(&i,2 - 0^,2) 

Ci,l(ai,2 - 6i,2) + Ci,2(6t,l - Oi,l) 
Ci,2(Oi,l - + Ci,l(h,2 - Oi,2). 

We did the calculation of the scalar integral in three steps: 

1. The ^-integration, followed by a shift in Zo : ^1 to have the Zi in a standard 
form. 

2. A further shift Zi — >• Zi = W + Iq linearizes the z^s and allows for an application 
of the residue-theorem to do the /o-integration. 

3. Interpretation of the remainig integration as a representation of a TZ function. 

Arbitrary tensor integrals involve additional powers Zq 1 /™ 2 /^ 3 in the numerator so 
that we have to modify the three steps as follows: 

1. The Zj_-integration now gives — B(l + 71,3 — 6, — n% + e)(— Zi)~ e+Tl3 . The shift 
in Iq and l\ now introduces a (double) binomial sum of powers from Z§Z ° to 
Zq 1 /" 2 in the numerator. 

2. The further shift Zi — * Zo + Zi changes the powers in the numerator to have 

values from ZqZi to?Q 1+ ™ 2 Zi 1 . The application of the residue theorem inserts 
the residue also in the numerator, so that the remaining l\ integration involves 

,~0 r n 1 +2n 2 , 

powers Li . . . Li m the numerator. 

3. The remaining integrals are of the form J °° dy y ^ uv + ^ v ^^ so that 

we will end up with results of the form ~ B(k + 1, 2e — 2n 3 — ^ — 

K-2e+2n 3 +k(e -n 3 ,e- n 3 , 1; -y + , -y_, s/r). 

Applying reduction formulas for 7?.-functions all tensor integrals are now ex- 
pressible in a standard set of well-known functions Integrals involving higher 
powers of propagators can again be obtained via differentiation formulas. 

2.3 The Singular Case 

In evaluating the three- and higher n-point functions one counters more than just 
the C/V-singularities. There are on-shell and infrared singularities as well. These 
divergences always appear as parallel-space divergences. Here they enter the calcu- 
lations in the form of endpoint singularities. So they can be systematically handled 
by applying Hadamard's finite part to them pj. 



A i>2 = A iA = 
A,2 — — A',4 = 



G'1,1 — 


—Ci,3 


Ci,2 = 


— Cj,4 


A,i = 


= A,3 


A,2 = 


= A,4 
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3 2-loop Functions 



In this section we start with the calculation of massive scalar 2-loop functions. We 
start with the 2-point case as a more detailed example and present the results for 
the two different topologies which appear in the 3-point case. 



3.1 The 2-point Function 

Let us first introduce some notation. The expression for the normalized two-loop 
two-point function is 



According to our method of separating parallel and orthogonal space variables 
we rewrite the integral as 

I (l 2 ) = -^jd% Jd%f[±-, (10) 

where the inverse propagators Vi are now defined by 

Vi = l\ n - - m\ + ig 

V 2 = (ho + qo) 2 - li± - m\ + ig 

V 3 = (ho + ho) 2 - l 2 ± - l%± - 2h±h±z -ml + ig 

Va = (ho - qo) 2 - l 2 ±_ - ml + ig 

V 5 = l\ a - l% ± -m\ + ig, 

and we choose to be in the rest frame of the exterior momentum q = q . 
The loop momenta h an d h satisfy 

h = (ho,h±) 

h = (ho,h±) 

h± := | h± | 

h± '■= I h± I 

h ■ h = h±h±z, 

where the space-like h± and h± span the orthogonal space. 

Note that we have to introduce an angular variable for the relative angle between 
h± and l 2 ±- We have to do so because we are not allowed to choose h± • h± = 0. 
This choice would fail to give the same results for DR and ordinary Riemann 
integration in the case of finite integrals. 

The integration over the orthogonal space involves six integrations. Three of 
them give the product of the volumes of the unit spheres S 1 and S 2 as a prefactor 
because only one non-trivial angular integration is involved in the scalar-product 
h ■ h- As a consequence we are left with only five non-trivial integrations, two 
(ho an d ho) from the parallel space and three (h±,h± and z) coming from the 
orthogonal space. We are left with the the following integrations to perform 

Q„2 POO />OC f' oo /*oo /*1 5 -. 

I(q 2 ) = ^^ dl w dho dli± dh±ljj 2 2± dz\\—. (11) 

71" J-oo J-oo Jo JO J-1 fj[ l J i 



G 



As usual g is small and strictly positive and the integrations are well defined 
for q 7^ 0. Because we do not Wick rotate we keep the i ^-prescription throughout 
the whole calculation. At the end of the calculation we will find the limit g — > to 
be well-behaved in every case. 

This 'art' of calculating complicated two- loop functions can be described as 
transforming the integrand so that it can be recognized as an integral representa- 
tion of special functions Q. Usual methods use Wick rotations, so an analytical 
continuation to physical values at the end of the calculation is always understood. 
As a consequence, to have sensible results, the final expressions must be given 
in terms of known special functions in order to manage the analytic continuation 
back. The zp-prescription is then re-installed in order to make this continuation 
unique. In our approach we avoid altogether the need to Wick rotate. No analytic 
continuation at the end of the calculation is necessary and sensible results can be 
obtained by numerical evaluation of the final integral representation. 

Let us proceed with the integrations (a detailed discussion can be found in [p"l|). 
To do so, we first integrate z, the cosine of the angle between l\ and I2. 

The integral becomes 



Expanding the logarithms and using symmetry properties of the integrand one 
finds with the help of the residue theorem the simple result 




x 



[log((/io + ho) 2 
log((Zi + ho) 2 



■ml+ig- (hx + hxf) 
mj+ig- (h± - h±) 2 )] ■ 




C 




1 

log(lOj + Wj + Wk) 

—m\ + m\ — q 2 — 2lioq 
—ml + m\+ q 2 — 2Z 2 o<3, 



or in dimensionless variables x = ho/q, y — ho/l 





x [£325 — £315 + £314 
4 




(12) 



log(wj +Wj+ Wk) 

-z\ + z\ - 1 - 2x 
-zl + zl + 1- 2y, 



where zf 



mf/q 2 replaces mf in each Wi'. 
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w 2 = 




+ i) 2 


-4- 


h ig 


w 3 = 




+ y)2 


-4- 


1- ig 


U>4 — 




-1)2 


~ z\ - 


- ig 


W 5 = 


\A 2 


-4- 


V ig. 





Note that for each Wi we have < arg(?«i) < ir/2 so that the limit g — > turns out 
to be well defined. A detailed discussion of this result can be found in This 
result was tested and confirmed by analytical and numerical checks and was widely 
used in comparison with other numerical routines J^, || and asymptotic expansions 
I- 

3.2 The 3-point Functions 

In the following we will present a integral representation for 2-loop 3-point functions 
similar to the one for the 2-loop 2-point function. 

In the 3-point case, we have two scalar master diagrams according to the two 
possible different topologies, Fig. (Jl|) and Fig.(||). We will find the amazing result 
that for both topologies a threefold integral representation can be given over the 
same simple functions as in the 2-point case. Only the number of terms is different 
(from 36 for Fig.(|) to 80 for Fig.(0)). 

For 3-point functions additional problems may arise from on-shell singularities. 
As we will see our method offers a very systematic handling for these divergences. 
Because only endpoint singularities arise in our calculations, it seems very ap- 
propriate to extract these singularities via principle value prescriptions. A brief 
introduction to this method is given in fl^| . 

It is an amazing consequence of our method that the crossed diagram Fig.dll) is 




Figure 1: The Cross Diagram, the labels denote the propagators 

of the same order of difficulty as the parallel diagram Fig.(||). This is a consequence 
of the partial fractionizing, which assigns, roughly speaking, twice the number of 
terms to the crossed diagram as to the planar one. We begin with a short discussion 
of the planar function. 



8 



Figure 2: The Planar Diagram, the labels denote the propagators 

The generic function is a function of nine variables (six masses and three kine- 
matical variables). First we give the explicit expressions for the propagators: 

Pi = (la + qio) 2 mf + in 

P% = (la - <72o) 2 - (h - <?2i) 2 - l\ - m\ + irj 

P3 = (la + h) 2 - (h + h) 2 - l 2 ± - k 2 ± - 2l ± k ± z-ml + ir] 

Pi = (k - qio) 2 -k\-k\-m\-\- ir) 

P5 = (ko + Q2o) 2 - Oi + 92i) 2 - k\ - ml + irj 

Pe = k — k\ — k ± — m 6 + irj. 

Again z denotes the angle between l± and kj_ in the orthogonal space. The expres- 
sion for the whole integral is then 

p + 00 poo n+l 11 1 

I = C dlodhdkodkj dl ± dk ± dz p p p p p p ■ (13) 

We first integrate z and then make the shifts Iq — > Iq + l\ and fco — > ko + ki 
which locates both the cuts in l\ and k\ of the square-root in the result of the 
z-integration in either the upper or lower halfplanes, depending on the sign of 
(la + ko). Furthermore, all propagators become linear in Iq (resp. fco). So we can 
do without any difficulties the l\ and k\ integrations next using the residue theorem 
(details of all these calculations are in |Q ). 

The explicit form of the remaining propagators is now 

Pi = (la + h + qio) 2 - l{ - l± ~ + irj 

Pi = (lo + h - <72o) 2 - (h - <?2i) 2 - l\ - m\ + irj 

Pi = (ko +ki — qio) 2 — ki — k± — m\ + irj 

P 5 = (k + ki + q 2 o) 2 - (ki + q 2 i) 2 - k\ - m\ + it] 

P 6 = (k + k^ 2 ~ k 2 ~ k 2 ± -ml+irj. 

Again we use a partial fraction to simplify the structure of the above product of 
propagators. Each of the above propagators involve either I or k variables, but not 
both. So we partial fractionize the I and k propagators separately: 

1 1 

PlP 2 PiP 5 P6 = Pi -Pi 



1 

P 



1 

Pi 



11 11 

L [Pi ~ Ps)(P4 -P 6 )P~4 + {P 5 - P 6 )(P 5 - Pa) Pn 
1 1 



(P 6 - P 5 )(P 6 - P 5 



(14) 



This decoupling of I and fc-variables reflects the nested loop structure of the 
planar topology of our graph. For the crossed ladder graph we will not have this 
decoupling of variables. 

Again all possible differences p — Pj are linear monomials in l\ or k\ with poles 
on the real axis (if all ?7's are assumed to be equal). 

All of these poles give in contributions to the corresponding integrations. But 
there is a further pole coming from the one linearized propagator which remains 
in each partial fraction summand. This propagator has a pole which is not located 
on the real axis. This pole contributes if it is in the interior of the contour of 
integration, and this in turn is determined by the sign of (Iq + kg). For example, 
if (Iq + fco) > 0, we have to close the contour in the upper halfplane for both the 
l\ and fci-integrations. Then the pole of P\ in l\ contributes for Iq < —qiQ. This 
gives a further splitting in the Zo-integr&tion. 

A list of all the resulting monomials, sufficient for the construction of the full 
result, together with the corresponding REDUCE program, will be published else- 
where |Q. It can be also obtained from the author on request. 

Having done the l\ and k\ integrations, one ends with expressions of the fol- 
lowing form for the l\ and k\ dependence, 



^,2^2 l/y/Poly(^,fcj) n ,, 
o dl ^(a l ll + b l )(a k kl+b k y (15) 

where Poly(^, k\) is a polynomial expression in k\ of second degree 

Poly(x,y) = a(x - x + (y))(x - x_(y)). (16) 

Let us do the /^-integration next. Expressing the integrand in terms of 1Z- functions 
and applying a Landen transformation || one can express the result in terms of 
logarithms and square roots |l^| 

K-i(l/2, 1/2, 1; z+, Z-, z ) = 7£_i(l, 1; u + , u_) = , (17) 

w+ — U- 

where u+, u- are given by 

U + = \Jzl~ Z (Z + + ZJ) + Z + Z- + Zq^/Z+Z- 



zl - z (z + + Z-) + z + z_ - Z 0y /Z+Z-. 



(18) 



where the z\ variables are polynomial expressions in the remaining integration 
variables. 

In the above mentioned programs the explicit expressions for Poly(Zj_, kj_) are 
generated. From them, one can easily determine the resulting w + and U- expres- 
sions. So we end with a representation as a threefold integral for the generic case 
of six different masses. In all cases the Iq (resp. fco)-axis is divided in at most two 
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pieces through the splitting of domains which means that we have four domains to 
investigate. We have six summands from the partial fraction and six (three from 
the k\ integral times two from the h integral) contributing residues in each domain. 
This gives 36 terms, each one having its own domain of contribution. 

The crossed function is regarded to be much more difficult because of its non- 
planar topology. But we will see in a moment that again, with partial fractions 
and shifting of momenta, we can express it in terms of the same function as for the 
planar one. Only the number of terms increases. The crucial observation is that 
we can do a partial fraction decomposition in the two propagators involving both 
loop momenta. Let us first list the propagators explicitely: 

Pi = (lo + k - qw) 2 - (h + kif - l 2 ± - kj_ - 2l±k±z - m? + ir) 

P 2 = {Iq + h + q 2a ) 2 - (7i +fci +<72i) 2 -l 2 j_-k 2 j_-2l ± k ± z-ml + ir) 

P 3 = {lo - qw) 2 -lj-l 2 ± -rnl + in 

Pa = {ka + <22o) 2 - (fa + q 2 i) 2 - k\ - m\ + irj 

P 5 = ll-lj-ll-ml + if] 

P& = k — fa — — m 6 + ii]. 

Then we have 

1 _ 1 
~P\~p2 ~ Pi -Pi 
Pi -Pi = -2{q w + q m )(l + k ) + 2q 21 (h + fci) - mj + m\. (19) 

In the difference Pi — P 2 the orthogonal space dependence cancels out completely. 

In each summand we can make a shift so as to have the same structure after the 
z-integration just as for the planar case. Then we can proceed exactly as for the 
planar diagram. The one crucial difference lies in the extra pole l/(Pi — P 2 ) ■ This is 
a pole which mixes I and k variables. Nevertheless, the li, ki and l\ integrations can 
be done in a similar manner, only the book keeping of the contributing expressions 
and the relevant splittings in the lo and ko become rather tedious. This book 
keeping is done again with the help of RED U CE and is included in |Q . 

It is remarkable that we end up with integral representations of the same kind 
as in the planar case while only the number of terms increases. Explicitely, we have 
three possible residues after the integration. Inserting the Pi — P 2 pole, we have 
four poles in the following ki integration. The other two poles in li only lead to 
three poles only by the following fci integration. So we end with 4 + 3 + 3 = 10 
contributing residues. For the crossed function, we have 2x2x2 = 8 summands 
from the partial fractions (instead of 2 x 3 = 6 in the planar case). Altogether we 
have 80 terms which contribute. 

Let us mention here that possible on-shell divergences which may appear will be 
handled by using Hadamard's finite part. A detailed demonstration of this method 
will be presented elsewhere jl5| (see also [[l3| [l2|, |Io|). 

Future work will be devoted to implement this Hadamard subtraction prescrip- 
tion in our programs. Note that for the singularities that have been extracted one 
fixes one or two integration variables, so that we have only one- or twofold integral 
representations for the coefficient of the singularity, which can be evaluated at least 
numerically. A list of the relevant coefficients in the on-shell case (with the ladder 
(resp. cross- ladder) propagators assumed massless) is already included in Q]. 
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4 Conclusions 



We have presented in examples methods which apply to all loop integrals which may 
appear in the Standard Model up to the 2-loop level. We hope that these methods 
add to the set of tools used in multiloop calculations. Apart from a remarkable 
systemization of 1-loop functions we have derived integral representations valid 
for planar and non-planar massive 2-loop 2- and 3-point functions. We calculated 
these results for the generic case, that is as functions of nine variables {qo,qi,qo ■ 
qi, six masses) in the three-point case. Accordingly, the results turn out to be 
complicated. Nevertheless they are found to have a very systematic structure 
which reflects itsels in a repeated pattern of typical terms. The extraction of 
possible infrared singularities is reduced to a book keeping of end-point singularities 
for both topologies. In future work, one must investigate how much CPU time is 
necessary for all the book keeping. We hope that the above integral representations 
turn out to be manageable with respect to these practical considerations. 

Also one can try to find easier integral representations if one treats simpler 
mass cases. There are enormous possibilities for the application of these results. 
Almost every calculation in the electro-weak sector involves integrals with different 
masses, but typically one can restrict ones attention to three different masses. So, 
if the above integral representation allows for an evaluation in a reasonable CPU 
time, this may open the door to a new era for calculations devoted to a test of the 
Standard Model at two-loop order. 
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A Hadamard's Finite Part 

We give a short summary of the properties of Hadamard 's finite part. Our presen- 
tation follows |l6| . We use this distribution to have a systematic regularization for 
on-shell and infrared singularities. 
Consider the integral 



where cj) belongs to an appropriate set of test-functions. Let fit") be locally inte- 
grable but singular in the neighborhood of t = a. 

Let us assume that it is possible to subtract from F(a, b, e) a finite linear combi- 
nation of negative powers of e and positive powers of log(e) such that the remaining 
function of e allows for a well-defined finite limit e — > 0+. Then, this limit is called 
Hadamard's finite part of /. It defines a distribution H(f). The functions / which 
allow for such a subtraction of the singularities are called pseudofunctions. 

To subtract the singularities in e explicitly, one uses a Taylor expansion of the 
test-function <p(t) = <j>{a) + (a)(t — a) + . . . at a. For a pseudofunction /, only a 




> 



(20) 
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finite number of terms is not integrable (in the limit e — ► 0+) when the expansion 
of cf> is inserted in Eq.(pO|). Let us assume that the first k terms in the Taylor 
expansion are ill-defined. Then the finite part is given by 



while the singular part is a divergent combination of negative powers of e and log(e) 
given by the integral of the first k terms. 

The only property of the test-function, that is used in the above definition 
of Hadamard's finite part is the fact that it is integrable in the interval (a, b). 
Therefore, we can define a regularization scheme by the use of the above procedure. 
It is applicable for every integrand which belongs to the class of pseudofunctions. 
This is sufficient for our applications because the wy-prescription guarantees every 
Feynman integrand to be a pseudofunction. 

The scheme assigns a linear combination of negative powers of e and positive 
powers of log(e) to the singular part. For the finite part, its result is Hadamard's 
finite part. The generalization to iterated integrals is straightforward, one mainly 
has to replace the above Taylor expansion in one variable by an expansion in several 
variables. 
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